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Abstract. It is an open question whether solutions of the Einstein- 
Euler equations are smooth enough to admit locally inertial coordinates 
at points of shock wave interaction, or whether “regularity singularities” 
can exist at such points. The term regularity singularity was proposed 
by the authors as a point in spacetime where the gravitational metric 
tensor is Lipschitz continuous (C 0,1 ), but no smoother, in any coor¬ 
dinate system of the C * 1 ' 1 atlas. An existence theory for shock wave 
solutions in C 0,1 admitting arbitrary interactions has been proven for 
the Einstein-Euler equations in spherically symmetric spacetimes, but 
C 1 ’ 1 is the requisite smoothness required for space-time to be locally 
flat. Thus the open problem of regularity singularities is the problem 
as to whether locally inertial coordinate systems exist at shock waves 
within the larger C 1 ’ 1 atlas. To clarify this open problem, we identify 
new “Coriolis type” effects in the geometry of C 0,1 shock wave metrics 
and prove they are essential in the sense that they can never be made to 
vanish within the atlas of smooth coordinate transformations, the atlas 
usually assumed in classical differential geometry. Thus the problem of 
existence of regularity singularities is equivalent to the question as to 
whether or not these Coriolis type effects are essentially non-removable 
and ‘real’, or merely coordinate effects that can be removed, (in analogy 
to classical Coriolis forces), by going to the less regular atlas of C 1 ’ 1 
transformations. If essentially non-removable, it would argue strongly 
for a ‘real’ new physical effect for General Relativity, providing a phys¬ 
ical context to the open problem of regularity singularities. 


1. Introduction 

A basic question for the shock wave theory in General Relativity (GR) 
is the regularity of the gravitational metric at shock waves. Shock waves 
are solutions of the Einstein-Euler equations in which the fluid density and 
velocity are discontinuous. Shock waves always form in solutions of the 
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compressible Euler equations whenever the flow is sufficiently compressive, 
urn- Shock waves introduce increase of entropy, time-irreversibility and 
loss of information into GR, and they also create discontinuities in the cur¬ 
vature tensors of space-time. Classical shock waves in non-relativistic gas 
dynamics are regularized by shock profiles when viscosity and heat conduc¬ 
tion are included, but the theory of dissipation is problematic in relativity 
due to the fact that parabolic equations introduce infinite speed of propa¬ 
gation, and modified theories of dissipation have been controversial as they 
typically are either not causal, or do not admit shock profiles (see [5] for 
references) Q Moreover, the discontinuities that appear in the zero dissipa¬ 
tion limit are replaced by steep gradients near the limit, and the essential 
issues of shock waves persist. Thus the Einstein-Euler equations have played 
a fundamental role in relativity both because they accurately model highly 
relativistic flows, and because, at a fundamental level, shock waves describe 
an accurate idealized limit which introduces dissipation into relativity with¬ 
out giving up causality. Even so, it remains an open problem as to the 
regularity of the gravitational metric at GR shock waves. Specifically, it is 
not known whether the space-time is smooth enough to admit locally in¬ 
ertial coordinates in which the metric is Minkowski and the derivatives of 
the metric vanish at points of shock wave interaction @ If space-time is not 
locally inertial in the zero dissipation limit, then the lack of existence of 
locally inertial frames in this limit would persist as an issue under pertur¬ 
bation by shock profiles because the jumps in first-order derivatives would 
get replaced by large second-order derivatives, and the irregularity of the 
metric would simply be spread out over a small region which would appear 
to be a regularity singularity in the far field limit. 

A basic existence theory for shock wave interactions in the Einstein-Euler 
equations based on Glimm’s method was established in [6j for spherically 
symmetric spacetimes, and interestingly, the methods are only sufficient to 
prove existence for gravitational metrics g which are only Lipschitz con¬ 
tinuous. It is not known whether the metrics associated with solutions in 
[S|, or any more general shock wave solutions, can be smoothed from C 0,1 
to C 1 2 ’ 1 by coordinate transformation at points of complicated shock wave 
interaction]! The metric regularity C 1,1 is the minimum regularity that 
guarantees a metric admits locally inertial coordinate frames, and for the 
weak and strong formulation of the Einstein equations to be equivalent, M- 


1 See [5] for a new theory of dissipation based on including relativistic viscosity and heat 
conduction parameters to obtain a symmetric hyperbolic regularization of the classical 
relativistic Navier-Stokes-Fourier equations for pure radiation, such that the resulting 
equations are causal and dissipative, and such that all shocks admit shock profiles. 

2 We say a point of shock wave interaction is a point where multiple waves interact in 
any complicated fashion, so long as one of them is a shock wave. 

'^Here C 0,1 denotes continuous with Holder derivative one, (i.e., Lipschitz continuous, 
0), so for metrics g € C 0 ’ 1 , first derivatives of the metric suffer a jump discontinuity at 
shocks, while for g £ C 1 ’ 1 the second derivative suffers a jump discontinuity at shocks. 
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Moreover, within the C 1)l atlas, the condition that solutions be free of delta 
function sources is a covariant condition, and this appears to be the weakest 
atlas with this property, appropriate for shock wave solutions in GR, c.f. 

B51- 

To set up a framework for addressing the question as to whether metrics 
associated with shock waves can be smoothed from C 0,1 to C 1,1 by coordi¬ 
nate transformation, we begin below by proving that if g is a shock wave 
solution which is only C 0,1 regular in one coordinate system, then the metric 
cannot be smoothed to C* 1,1 within the atlas of smooth (say, C 2,1 ) coordinate 
transformations, the atlas usually assumed in classical differential geometry. 
Therefore such metrics do not admit locally inertial coordinate frames within 
the smooth atlas, (c.f. Theorem l3.2l belowl. However, the Einstein equations 
remain consistent in the weak sense when the smooth atlas is extended to 
the larger atlas of C 1,1 coordinate transformations, [IB]- The Jacobians of 
C 1,1 transformations are only C 0,1 , and thus C 1,1 transformations have the 
required properties to potentially eliminate jumps in the metric derivatives 
at shocks. Thus the more singular atlas of C 1,1 coordinate transformations 
is the atlas that holds the possibility of lifting the metric regularity to C 1,1 , 
mm • Thus a most natural open question is whether the gravitation metric 
can be smoothed from C 0,1 to C 1,1 by C 1,1 coordinate transformations, and 
this addresses the locally flat character of space-time at GR shock waves. 

Following our work in |T2|, we define a regularity singularity as a point in 
space-time where the metric does not admit locally inertial frames within the 
C 1 ’ 1 atlas, and it is an open problem as to whether regularity singularities 
can be created by shock wave interactions in GR. The authors know of no 
physical or mathematical principle that can rule out regularity singularities 
in the Einstein-Euler equations ahead of time, and only mathematical proof 
can ultimately resolve the issue. 

The starting point for our analysis here is the celebrated paper m by 
Israel, in which he shows that at smooth co-dimension one shock surfaces in 
n dimensions, the gravitational metric can always be smoothed from C 0,1 
to C 1,1 by introducing Gaussian normal coordinates (GNC) at the shock. 
The transformation to GNC is a C 1,1 transformation [16] . GNC are only 
defined for single, non-interacting shock surfaces and do not exist for the 
more complicated C 0,1 metrics constructed in the Groah-Temple framework 
[6]. Since Israel’s result in (8), it has been unknown whether the regularity 
of the gravitational metric can always be smoothed from C 0,1 to C 1,1 by 
coordinate transformation at shock wave interactions. Thus it remains an 
outstanding open problem as to whether space-time is always locally flat at 
points of shock wave interaction in GR. 

The first extension of Israel’s theorem to the more complicated setting 
of shock wave interactions was accomplished by the authors in their recent 
papers mm • The proof demonstrates in the first case ever that the grav¬ 
itational metric can always be smoothed from C 0,1 to C 1,1 at a point of 


4 


M. REINTJES AND B. TEMPLE 


interacting shock waves in GR, namely the case of regular shock wave inter¬ 
action in spherical symmetry between shocks from different characteristic 
families. The proof introduces an explicit physical procedure for finding the 
coordinates that display the physics in their simplest form, replacing the 
GNC construction in an essential new way, based on solving a non-local 
hyperbolic-type system of equations. But the argument is tailored to the 
specific case of two interacting shock waves, and there are many miracles 
in the proof, places where the Rankine-Hugoniot jump conditions come in 
to surprisingly make an apparently over-determined system barely solvable. 
Thus it is not clear whether or how this proof can be extended to more 
complicated shock wave interactions. For complicated shock interactions, 
including the more complicated C 0,1 solutions that exist by [6], or for com¬ 
plicated asymmetric shock interactions in (3 + l)-dimensions, the question as 
to the locally flat nature of space-time and whether regularity singularities 
can be created by shock wave interactions, remains an open problem. 

Our work on the problem has led us to conjecture that if regularity singu¬ 
larities actually exist, the structure of space-time at the singularity would be 
essentially determined by its structure within the smooth atlas alone. That 
is, at a regularity singularity, the larger C 1,1 atlas would offer no essential 
improvement in the regularity of the gravitational metric over and above 
that observed in the C 2,1 atlas. Thus to begin studying the implications of 
regularity singularities, should they actually exist, we first determine prop¬ 
erties of a Lipschitz continuous space-time metric within the smooth C 2,1 
atlas in a neighborhood of a point on a single shock surface. These proper¬ 
ties extend easily to the case of shock wave interaction. Our purpose in this 
paper, then, is to establish physical implications of the assumption that the 
gravitational metric is no more regular than its regularity within the smooth 
atlas. 

In Theorem 13.21 we prove that, restricting to the smooth atlas for a met¬ 
ric Lipschitz continuous across shock surfaces, the closest one can get to a 
locally inertial coordinate system is one which is approximate locally inertial 
in a natural sense we make precise. We then characterize what we interpret 
as Coriolis type effects in approximate locally inertial coordinates, effects 
which arise from terms which only vanish in a true locally inertial coordinate 
system, should one exist. These Coriolis type effects are analogous to clas¬ 
sical Newtonian Coriolis forces which are due to terms in the gravitational 
force law which arise from the rotation of the earth, but would vanish in a 
true locally inertial coordinate system. (Keep in mind that classical Corio¬ 
lis forces would be treated as real until inertial coordinates, which remove 
them, are identified.) In Section [5l we identify the Coriolis terms in the ge¬ 
odesic equations in approximate locally inertial frames. In Theorem 16.11 of 
Section O we derive a canonical form of the linearized Einstein equations in 
approximate locally inertial coordinate systems and use this to identify the 
Coriolis terms associated with gravity waves. In Theorem 16.31 we use our 
formulation of the linearized Einstein equations to prove that these Coriolis 
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terms are nonzero and cannot be removed by coordinate transformation to 
any approximate locally inertial frame. This quantifies the contributions to 
the scattering of gravitational radiation in approximate locally inertial coor¬ 
dinates. The main results are summarized in Theorem 15.11 Theorem 15.11 is 
interesting in its own right, because it describes how far a metric Lipschitz 
continuous across shock surfaces is from being locally inertial within the 
smooth atlas, in terms of gravitational radiation. Within this context, the 
open problem of regularity singularities, then, is the problem as to whether 
locally inertial coordinate systems and the essential regularity of space-time 
can be improved upon by extending the smooth atlas to the larger C 1,1 
atlas. We conclude that, if no such improvement exists, then the quantifi¬ 
able effects produced by these non-removable Coriolis terms are physical 
implications of regularity singularities. Q 


2. Preliminaries 

In General Relativity, the gravitational field is described by a Lorentzian 
metric g of signature (—1,1,1,1) on a four-dimensional spacetime manifold 
M. We call M a C^-manifold if it is endowed with a C^-atlas, a collection 
of four-dimensional local diffeomorphisms that are C k regular from M to 
R 4 . A composition of two local diffeomorphisms x and y of the form xo y -1 
is referred to as a coordinate transformation. 

Our index notation for tensors here sometimes uses indices to determine 
the coordinate system, e.g., T„ denotes a (1, l)-tensor in coordinates x fl 
and Tp denotes the same tensor in coordinates x a . We use the Einstein 
summation convention whereby repeated up-down indices are summed over 
all values for the given indices. Tensors transform by contraction with the 
Jacobian and the inverse Jacobian, which we denote by 

J“ whenever there is no confusion. In particular, the metric transforms as 
9fiu = J^Jvda/3- Tensor-indices are raised and lowered with the metric g /lv 
and its inverse g^ v . 

The time evolution of a gravitational field in general relativity is governed 
by the Einstein equations [3] 

G ij = nT ij , (2.1) 


4 The assumption that the curvature tensors contain no delta function sources is crucial 
for the conjecture as to whether regularity singularities exist in GR. Shock surfaces, by 
definition, are weak solutions which contain no delta function sources. Interfaces which 
contain delta function sources were introduced in Israel’s theory of thin shells, and these 
include the domain wall between the false vacuum and the true vacuum in Guth’s cele¬ 
brated theory of inflation, ]T|. Metrics with interfaces containing delta function sources 
lie in C 0 ’ 1 , but cannot be smoothed to C 1 ’ 1 within the C 1 ’ 1 atlas at the start because C 1 ’ 1 
metrics have classical curvature tensors which contradict the existence of delta function 
sources. But, when the curvature tensor is free of delta function sources, there appears to 
be no physical reason for this loss of regularity. 
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a system of 10 second-order partial differential equations that relate the 
metric tensor gij to the undifferentiated sources through the Einstein 
curvature tensor 

G ij = R ij - ~Rg ij , (2.2) 

a tensor involving second derivatives of g. The Ricci tensor, Rji, is the trace 
of the Riemann tensor, and the scalar curvature, R, is the trace of the Ricci 
tensor, that is, 

Rji = g lk R ijk i and R = g jl R j i. 

The Riemann tensor is given in terms of the Christoffel symbols T] ■ by the 
formula 

Rijkl — T ijk,l T'ijljk T g ^ (r ial^pjk pjl') j (2.3) 

where 

Fuj = - (.guj + gjij - g ij} i ) and Y k ^ = g kl Tnj . (2.4) 

and k = ~^G is the coupling constant which incorporates Newton’s gravita¬ 
tional constant G and the speed of light c@ Here is the energy-momentum 
tensor. For our methods we do not need to specify the matter sources, but 
our motivation comes from the case when is the energy-momentum ten¬ 
sor for a perfect fluid, 

T tj = (p + p)v,iUj + pg t j , (2.5) 

where p is the energy density, U{ the unit 4-velocity, u l iii = — 1, and p the 
pressure. 

Conservation of energy and momentum enter the Einstein equations through 

T ij ;j = 0 , (2.6) 

which reduces to the relativistic compressible Euler equations in flat space- 
time, and follows from the divergence-free property of the Einstein equa¬ 
tions, G lJ ^ = 0, a property built into G at the start as an identity following 
from the Bianchi identities of geometry, [20]. Here as usual, semicolon de¬ 
notes covariant differentiation 


= ^.+ny, 


where T] ■ denote the Christoffel symbols associated with metric g, defined 
in (12.41) . In the special case of a perfect fluid, equations (I2.ip and (12.611 form 
the coupled Einstein-Euler equations; a system of second-order differential 


'’’Note that because the curvature tensor is anti-symmetric in ( k, l) it requires a choice 
of sign, and this choice is not uniform in the literature. Here we use Weinberg’s convention 
120) . also used by the tensor package of MAPLE, and with this sign convention for the 
curvature tensor, the gravitational constant k is negative, i.e., k = — The careful 
reader should beware that Hawking and Ellis, as well as the MAPLE geometry package, 
use the opposite sign convention, defining R to be minus our R , resulting in k = 

Note that Weinberg’s convention 77 = (1, — 1,— 1,— 1), different from ours, has no bearing 
on this choice of sign. We also point out that our methods apply essentially unchanged 
for non-zero cosmological constant. 
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equations for the unknown metric Qij . coupled to the fluid variables p, p and 
v?, a system which closes upon specification of the equation of state, [20] • 
In special relativity the spacetime metric is taken to be g tJ = rjij where 
r]ij = diag(— 1,1,1,1) is the Minkowski metric, in which case (12.61) reduces 
to the relativistic compressible Euler equations, a system of conservation 
laws in which it is well-known that shock waves form out of smooth initial 
data whenever the flow is sufficiently compressive. Shock waves are discon¬ 
tinuous solutions that solve the Euler equations weakly, in a distributional 
sense, mm- Across a smooth shock surface £, the Rankine-Hugoniot jump 
conditions hold, 

[T ij ]nj = 0, (2.7) 

where [/] = fh ~ fn denotes the jump from right to left in function / across 
E, and rij is the surface normal. In particular, for smooth shock surfaces, the 
jump conditions (12.71) are equivalent to the shock wave solution satisfying 
the weak formulation of (12.6|) across X, c.f. m- 


3. The smooth (C 2,1 ) atlas at points of shock wave interaction 

AND APPROXIMATE LOCALLY INERTIAL FRAMES 

At a point p where the space-time metric is at least C 1,1 , the Riemann 
normal coordinate construction implies that one can always transform the 
metric to locally inertial coordinates at p , coordinates for which gij(p ) = 
rjij = diag{— 1,1,1,1) and gap^ip) = 0 so that the metric can be Taylor- 
expanded and written as 

9 a p = Va/3 + 0(5 2 ), (3.1) 

where <5 is the 4-dimensional (Euclidean, non-covariant) coordinate distance 
to p. That is, 


5(q) = 5(x(q) -x{p)) 


\ 




Ct=0 



(3.2) 


In Theorem 13.21 below we show that the closest you can get to a locally in¬ 
ertial coordinate frame within the atlas of smooth (at least C 2,1 ) coordinate 
transformations for a Lipschitz metric at a point of shock wave interaction 
is an approximate locally inertial coordinate system, a term we make precise 
in the following definition. 


Definition 3.1. We call a coordinate system x a an “approximate locally 
inertial coordinate system” at a point p if the metric takes the form 


9 = V + 9 

in a neighborhood V of p where g is a Lipschitz continuous symmetric tensor 
vanishing at p, satisfying the condition that there exists an open set U C 
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V containing p in its closure, constants M,M > 0, and indices a,/3 ,7 G 
{0,3}, such that g is smooth in U, and 

\9a/3,y \ > M in U, (3.3) 

and, as a consequence of the Lipschitz continuity ofg, 

|<M(«)| < MS(q) inU, (3.4) 

where 5 is defined in h3.$\) . 

Condition (|3.3|> says that the shock wave interaction at p is sufficiently 
localized that there still exists an open set U with p in its closure in which the 
metric is smooth, but the derivative of the metric has a lower bound within 
this set, because the metric derivative takes a minimum jump on some shock 
wave in the neighborhood. This is a very weak structural condition which 
should easily be met at a point of finite shock wave interaction, including 
the shock waves which have been simulated in m, by choosing U as one of 
the regions between two adjacent shock curves intersecting in p. The metric 
smoothness in U makes our analysis in Section [ 6 ] feasible. 

Our next theorem, gives a precise sense in which approximate locally in¬ 
ertial frames are the closest you can get to an actual locally inertial frame 
within the smooth atlas in a neighborhood of a point on a smooth shock sur¬ 
face across which the metric is only Lipschitz continuous. The result extends 
from single shock surfaces to points of multiple shock wave interaction. 

Theorem 3.2. Assume a gravitational metric g^ v is Lipschitz continuous 
across a smooth co-dimension one (shock) surface E in a given coordinate 
system in the sense that the metric is smooth (at least C 1,1 ) away from 
E and smooth tangential to E, and such that there exists a constant M > 0 
and indices such that the jump in the derivative, \g flv ,a], across the 

shock surface satisfies 

\\g^,a]\>M. (3.5) 

Then the following holds: (i) For any coordinate system x a that can be 
reached within the C 2,1 atlas, the transformed metric g a p = JaJJjg^u is 
Lipschitz continuous but no smoother and there exist indices a, (3 ,7 such 
that across the shock surface 

|[<M, 7 ](9)| >||J- 1 («)||- 3 M, (3.6) 

for all q G E inside the coordinate patch and where || • || denotes the operator 
norm induced by the maximum norm on R 4 , (c.f. (13.81) below), (ii) Given 
a point p G E, you can always find a coordinate transformation within the 
C 2,1 atlas such that g a p(p ) = p a p, and every such coordinate system is an 
approximate locally inertial coordinate system in the sense of Definition ^. 1\ 
but never exactly locally inertial. 

®In all of this paper g G C 1:1 (U ) suffices, except for Section l6.ll where g G C 2 ’\U) 
suffices. 
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The condition (13.61) shows that M is a uniform bound over all coordinate 
transformations with bounded Jacobian. Thus assuming the shock strength 
is on the order of the jump in the derivatives in the original coordinates, M 
is then on the order of the shock strength, giving it an invariant physical 
meaning. 


Proof. To prove (i), we follow the idea leading to the smoothing condition, 
first introduced in [T 2 ], which lies at the heart of the method in [111 13]: The 
covariant transformation law of the metric is given by g a p = JaJ^g^v, where 
Ja = ffs- is the Jacobian of the coordinate transformation. Differentiating 
this transformation law with respect to a (we assume the index a belongs to 
the coordinates x ti ) and taking the jump of the resulting expression across 
E leads to 


[9a/3,c r] — {[Ja,a\J(3 + Jq [Jp,<r\)9/J.u + Jtf Jp \9y.v,a\, 

for all indices a, [5 and a. The above equation holds point-wise on E and 
the first term on the right hand side vanishes since the Jacobians are con¬ 
tinuously differentiable. Thus, contracting the resulting equation with the 
inverse Jacobians J“ and J„ gives 

jP?J2l9a^\ = [<W ]■ 

Taking the maximum of the absolute value of all components of the previous 
equation, we get 

max | J“ JP Jf[g a ) 9 7 ] | = max | \g^ a \ \ > M, (3.7) 

/z,^,cr 1 ^ 1 p, v,a 1 1 


where the lower bound of the right hand side follows from (13.51) . The left 
hand side can now be bounded from above as follows: Let || • || denote the 
induced operator norm (taken point-wise for (/ £ E), that is, for the inverse 
Jacobian, (J -1 )^ = J“, we have 


|| J 1 ( 9 )|| = sup \ max | Jn(q)w a \ w a £ M, a € {0,..., 3} and max \w a 

l f-i M a 

from which, for fixed n,u, a, we successively get 

< ||J _1 || • maxlJ^JJ[ 9q/3>7 ]| 

< ||«/- 1 || 2 -max|jJ[ 5a/3 , 7 ]| 

a,p 

< ||«/" 1 || 3 -max|[ ffQ/3 , 7 ]|. 

ot,p ,7 


-*}■ 

(3.8) 


Using the above inequality to replace the left hand side in (13.71) . we obtain 

max | [g a p 7 ] | > ||j - 1 || _ 3 M, 

which proves (13.61) for some indices a, /?, 7 , completing the proof of (i). 

We now prove (ii). By the symmetry of g^ u ( p ) there exist an orthonormal 
basis of eigenvectors at p. Choosing the Jacobian on a neighborhood around 
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p to be the constant matrix that maps the coordinate basis of x ,J ' to this 
orthonormal basis at p, we achieve a coordinate system x a within the smooth 
atlas for which g a p(p ) = r) a p. We now prove that x a is an approximate 
locally inertial frame. Defining g = g — 77 , the Lipschitz continuity implies 
(|3.4|) with M being the Lipschitz constant of g at p. The inequality (I .'1.3 1) 
follows from the metric smoothness away from the shock surface together 
with (13.61) . where we can take U to be the open set to the left or right of £. 
This completes the proof of the Theorem 13.21 □ 


4. Regularity Singularities 

Assume a Lipschitz continuous gravitational metric whose Riemann cur¬ 
vature tensor is bounded with discontinuities containing no delta function 
sources, c.f. D2I- 

Definition 4.1. We say that a point p in space-time is a regularity singu¬ 
larity if there does not exist a locally inertial coordinate system at p within 
the C 1,1 atlas. We say a point p is a weak regularity singularity if there 
exists locally inertial frames at p, but the metric is not C 1,1 regular in a 
neighborhood of p, in any coordinate system of the C 1,1 atlas. 

We conjecture that if a regularity singularity exists, then the C 1,1 atlas 
would offer no improvement to the smooth atlas regarding locally inertial 
frames, ISCED US]. Theorem 13.21 implies that regularity singularities ex¬ 
ist when the atlas is restricted to the class of smooth coordinate transfor¬ 
mations. The open question, then, is whether weak or strong regularity 
singularities exist with respect to the C 1,1 atlas, at points of shock wave in¬ 
teraction. Note that if a metric can be smoothed to C 1,1 in a neighborhood 
of p, it must then admit locally inertial frames at p, but it is an open prob¬ 
lem as to whether a metric might admit locally inertial frames at p, but not 
be smoothable to C 1,1 in a neighborhood of p. Thus we make the distinction 
between regularity singularities and weak regularity singularities. 


5. Coriolis Type Effects in Approximate Locally Inertial 
Coordinates and Statement of the Main Result 


In a locally inertial coordinate system centered at a point p where the 
spacetime is locally flat, the motion of test particles in free-fall will follow 
geodesics which are straight lines to errors quadratic in coordinate distance 
to p. because the Christoffel symbols vanish at p, i.e., T\- (p) = 0. Thus, 
neglecting these second-order errors, the geodesic equation is, 


d 2 x l 

It 2 


= V 1 ; ,i J = 0, 


where t parameterizes the curve, c.f. m- In contrast, if p is a regularity 
singularity, then by (13.31) . in each approximate locally inertial coordinate 
frame, there exists components of the metric whose derivatives are bounded 
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from below by a constant M > 0 in an open set U containing p in its closure. 
Since we can always solve for the derivatives of the metric in terms of the 
Christoffel symbols at a point via the formula, 

5o/3,7 9ap,"/ l^a /?7 T I)fl7cn 

there exists some set of indices, which for simplicity we again label as a, (3, 7 , 
such that 

Kp\ > M /2, (5.1) 

in the open set U. To avoid technicalities, (possible cancellations), assume 
for simplicity that a = 0 = /3. Then, choosing initial data for the geodesic 
such that x° = 1 and x l = 0 for i = 1,3, the corresponding geodesic 
equation for the fixed indices a = 0 = /3 and 7 becomes initially 


dr 2 00 ’ 

thereby isolating the presence of an acceleration in the geodesic on the order 
of |rgo| > M/2 throughout an open set U in each approximate locally iner¬ 
tial coordinate frame. We refer to the non-vanishing terms T^ as Coriolis 
terms. Since by Israel’s theorem points on single shock surfaces admit lo¬ 
cally inertial coordinates, such Coriolis terms and the resulting acceleration 
effects can be removed at points p on single shock surfaces, but by definition 
such effects could not be removed if p were a regularity singularity. 

Thus, at a regularity singularity, the non-removable Coriolis terms create 
geodesic deflection and scattering effects which (by definition) could not be 
removed by coordinate transformations. That is, in case spacetime is locally 
inertial, the trajectory of a particle following a geodesic curve can locally be 
approximated by a straight line with second-order error terms (in coordinate 
distance), which are determined by the curvature of spacetime. In case a 
regularity singularity is present at p, (larger) first order errors, which are due 
to the Coriolis terms, deflect the possible motion of the trajectory beyond 
what would be expected by the second-order errors (and hence curvature) 
alone. Since the Coriolis terms are of order M , the order of the shock 
strengths, they would represent physical effects at a regularity singularity 
should one exist. 

The purpose of the remainder of this paper is to characterize the non¬ 
removable Coriolis acceleration terms that appear in the linearized Einstein 
equations associated with a given approximate locally inertial coordinate 
system in the sense of Definition 13.11 In particular, shock wave interactions 
will typically form in dense matter where geodesic motion of test particles 
would be highly obscured by the matter, but, gravity waves are weakly 
interacting with matter, and the scattering of such waves at regularity sin¬ 
gularities could be an effect that would in principle be measurable. In the 
next section we derive a canonical form for the linearized Einstein equations 
for gravity waves in approximate locally inertial coordinates, accomplished 
in Theorem 16.11 The main technical problem is to identify non-zero terms 
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first-order in the derivatives of g , which survive in the linearized equations 
after all cancellation is accounted for. Indeed, since the curvature tensor has 
no delta function sources, we know that delta function sources in the second 
derivatives must all cancel out in the curvature tensor in an approximate 
locally inertial coordinate frame. This then begs the question as to whether 
there is a corresponding cancellation in the gravity wave equations that 
makes it so that nonzero first-order derivatives have no physical observable 
effect due to cancellation as well. We answer this question in Theorem l6.3l bv 
proving that no such cancellation occurs. Combining Theorems 16. II and 16.31 
we obtain our main theorem which states that no such cancellation occurs, 
and the linearized Einstein equations cannot be reduced to the Minkowski 
wave equation at p in any approximate locally inertial frame. (We always 
assume a background metric with curvature tensor bounded and free of delta 
functions sources.) 

Theorem 5.1. For each approximate locally inertial coordinate system cen¬ 
tered at a point p, and for each perturbation of the gravitational metric of 
the form r/ + g + eh, there exist a wave gauge such that, to leading order in 
e > 0, the linearized Einstein equations take the form 

^□r jhji = K [fjl ~ ■ ll.ilif + 

(5.2) 

a2 a 

where \3 V = r/ pt7 q xP q x <t is the wave operator, T is the matter source and 
R is the bounded curvature for the background metric g = rj + g, T is the 
contribution to the matter sources due to the perturbation h, and the terms 
Cij are the new Coriolis terms which involve first derivatives of g and hence 
would vanish in an actual locally inertial frame, were one to exist. Moreover, 
for each locally inertial frame, equations \5. 2\) always admit solutions h and 
indicies j and l such that 

\Cji(h)\ > -M (5.3) 

in some open set containing p in its closure 0 

Since a gravity wave at a locally inertial point of spacetime evolves to 
leading order by the pure wave equation of Minkowski spacetime, the term 
C(h) supplies the accelerations which characterize the scattering of gravity 
waves by the singularity. Since C(h) depends on the wave h as well as 
on the choice of coordinate system, the challenge is to prove that in every 
approximate locally inertial frame, there exists an h where the magnitude of 
C is on the order of h, and this is the substance of the Theorem. It establishes 
that regularity singularities would create Coriolis type accelerations that can 
not be removed by coordinate transformation due to the non-existence of 

~That is, the regularity singularity creates non-removable accelerations in each approx¬ 
imate locally inertial coordinate system within the class of gravity waves that meet the 
wave gauge in that coordinate system, a physically verifiable condition. 
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locally inertial frames. We thus propose that if regularity singularities exist, 
this is a physical scattering effect of the gravitational held due essentially to 
the lack of regularity in the underlying spacetime geometry. 

6 . Linearizing the Einstein Equations at a regularity 

SINGULARITY 

In this section we formally derive the linearized Einstein equations in 
approximate locally inertial coordinate systems as defined in Dehnition l3.ll 
The main issue here is to incorporate the hrst-order derivatives of the metric 
into the linearized equations, (c.f. m for a detailed description of the 
linearization procedure around the Minkowski metric). 

To start, assume an approximate locally inertial coordinate system x 3 at 
p where 

9ij = Vij + 9iji (6-1) 

where g is at best Lipschitz continuous, and we always assume the curva¬ 
ture tensor is bounded and free of delta function sources. Gravitational 
waves propagating through this background spacetime are represented by 
a symmetric tensor hjj such that the perturbed spacetime metric takes the 
form 

9ij — 9ij T 9ij T c hjj , (6.2) 

for some sufficiently small constant e > 0. In the next paragraph we in¬ 
troduce the wave gauge, and linearize the equations in this gauge. This is 
accomplished formally by substituting (|6.2I) into the Einstein equations in 
the wave gauge, assuming the neighborhood size 5 is order e, and discarding 
terms of order e 2 . The linearized equations (|5.2D then emerge as the equa¬ 
tions, linear in h, which result at order e. We then prove that not all Cji can 
vanish in a given locally inertial coordinate system, and since the equations 
are linear in h, we conclude that there exist solutions h which experience 
nonzero accelerations due to the Coriolis terms Cji- In this argument we do 
not address the issue of the validity of neglecting the 0(e 2 ) terms based on 
the smoothness of the perturbations h. For example, in the case of two inter¬ 
acting shocks, a gravity wave perturbation eh would create a small change 
in the positions of the shock waves encoded in the metric g , and in such a 
case the derivatives of the perturbation eh would not necessarily be small. 
However, as long as we can restrict to a subset of the original neighborhood 
in which the perturbation h is smooth, the linearization procedure would 
be valid in this neighborhood, and the effect of the Coriolis terms on h 
would be observed in that neighborhood. Thus, for our purposes here, we 
are content with applying the formal linearization procedure and not con¬ 
cerning ourselves with the actual regularity of the perturbations h. In line 
with this, we here understand second order derivatives of the background 
metric tensor g l3 , which is C 1,1 regular in the region under consideration, in 
a point-wise almost everywhere sense. 


14 


M. REINTJES AND B. TEMPLE 


6.1. The Wave Gauge. To linearize the Einstein equations in approximate 
locally inertial coordinates, we introduce the wave gauge condition, (also 
called the harmonic gauge condition, MU!!), on solutions of the Einstein 
equations. The wave gauge condition removes the gauge freedom from the 
Einstein equations by reducing the leading order terms to the wave operator 
(i.e., the D’Alambertian), and provides a canonical form for the Coriolis 
terms in the context of gravitational radiation. Differently from the standard 
linearization procedure, we start by introducing a wave gauge condition 
on the initial data for h which is propagated exactly by the full nonlinear 
Einstein equations, rather than by the linearized equation. In this gauge, the 
Einstein equations in the highest order derivatives produce the hyperbolic 
wave operator (r/* J + g %J ) ^ 7-577 of the background spacetime. Now, g is 
0(5 ) and 5 = 0(e) so only survives in the leading order derivatives 

at order 0(e) in the linearization procedure. The main point then is that 
derivatives of g are 0 ( 1 ) in approximate locally inertial coordinates, so they 
appear at order O(e) on lower order derivatives of the equations in wave 
gauge, the result being that the first-order derivatives of g only affect the 
lower order derivatives of the resulting linearized equations at order 0(e). 
These terms then constitute the sought after Coriolis terms. In particular, 
by this choice of wave gauge, it follows that we can only expect the wave 
gauge condition to be maintained within 0(e 2 ) errors under evolution by 
the linearized equations. 

Following the development in [2], we define the wave gauge condition by 

LOi = 0, (6.3) 


where 


for 


so that 


— £9 


!h: ij — ^ (hik,j T hjk,i h'ij.k ) , 


fT/. jj — Tfc ij E/. ij , 


(6.4) 


with Tpij denoting the Christoffel symbols of g = g + eh and Tkij the 
Christoffel symbols of the background metric g = 77 + g. As was shown in 
Choquet-Bruhat’s pioneering work on the existence theory of the Einstein 
equations, (16.31) is a condition on the Cauchy data only 12 El- That is, 
Choquet-Bruhat proved that whenever the initial data of a solution of the 
Einstein equations satisfies (|6.3I) . the Bianchi identities ensure that (|6.3I) 
holds in the whole Cauchy development. We explain this remarkable prop¬ 
erty of the Einstein equations in the following formal exposition. 

In her celebrated existence theory, to overcome the problem of the Ein¬ 
stein equations being degenerate hyperbolic, due to the constraint equations, 
Choquet-Bruhat’s strategy was to introduce the reduced Einstein equations , 
a modified hyperbolic version of the Einstein equations to which the Leray- 
existence theory can be applied directly. The idea is that the solution of 
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the reduced Einstein equations then solves the Einstein equations as long as 
it satisfies the wave gauge. For our purposes we take the reduced Einstein 
equations to be 

Hij\g\ = k ^Tij - ^g ij T ap g ,Tp \ , (6.5) 

DivgT = 0, (6-6) 

where in a given coordinate system x J the “reduced” Ricci tensor is defined 
in terms of the Ricci tensor and the gauge fields by 


Hij [g] — Rij [ 5 ] T tO i : j T !j 0 j.i , 


(6.7) 


and (|6.6I) expresses conservation of energy for the matter fields, which no 
longer is a consequence of the coupling to the Einstein tensor, c. f. mm- 

The system (I6.5D - (|6.6I) is a system for the pair (g t j , T tJ ), giving 14 equa¬ 
tions in 14 unknowns. To compare, the original Einstein equations consist of 
four constraint equations and six evolutionary equations for g t j. Assuming 
the constraint equations are satisfied initially, and assuming the six evolu¬ 
tionary equations hold, one can use the Bianchi identities to replace the 
constraint equation by Div g T = 0, resulting in ten evolutionary equations 
in 14 unknowns. This leaves four degrees of freedom upon which to impose 
gauge conditions. The idea of the reduced Einstein equations is to give up 
the gauge freedom. Namely, fix any coordinate system x 1 , then 1)6.511 and 
(|6.6I) give 14 equations in 14 unknowns in x l . The ten equations for the 
metric are second-order hyperbolic and the four for conservation are first- 
order hyperbolic, and there is a unique evolution of initial data. We show 
in the following that this evolution agrees with the evolution of the Einstein 
equations whenever the gauge condition (16.311 holds. 

To start, we first show that H\g\ is a hyperbolic operator. For this, it 
suffices to consider the terms of the Ricci tensor containing second-order 
derivatives. By (12.31) - (12.41) . these terms are given by 


9 ^ij[k,l\ — r)9 ( 9ik,jl 9jk,il 9il,jk + 9jl,ik ) 


Ak 


( 6 . 8 ) 


Ak , 



1 

1 Ak 

1 

1 9 

9jk,il 2 9ik,jl 

+ 9 

9il,jk 2 9ik,jl 


= 1 ( ra, 


9jk,i 2 9ik,j 


+ro j 


9il,k ^9ik t l 


where ^ g 9ji = g lk 9ji,ik denotes the wave equation based on g and 


^-ij[k,l\ — ^ijk,l r ijl,k■ 

A straight-forward computation starting from the definition of the Christof- 
fel symbols (12.41) shows that 

9 lk (vjk.i ~ 7>9ik,j) = 9 ap ^ja P = 9 ap (t'j'xp + efjo-p) > 
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where we used that I\,^ = f ij + e f ^. We now conclude by (16.41) that 

ddj = 9 ^ 9jk,i -^9ik,j\ 9 ' f'jifc. (6-9) 

Substituting (16.91) into (16.81) gives the leading order part of the Ricci tensor 
as 

— (wj,z + wzj) + l.o.t. 

= ^Ughji - (u jt i + Uij) + l.o.t., (6.10) 

where l.o.t. denotes sums of all terms not relevant to the hyperbolic structure 
of the Einstein equations for h, that is, sums of terms containing h and 
first derivatives of h, together with terms containing up to second-order 
derivatives of the background metric g. From (16.101) . it follows that 

Rij —\Zlghji (-jjjj T cu^j) 4“ l.o.t. 

and hence 

R-ij[ 9 ] — 4” l.o.t .. 

The conclusion then is that the reduced Einstein equations are hyperbolic 
and reduce to the Einstein equations whenever the gauge condition (16.31) 
holds, c.f. g]. 

In this context, we now outline Choquet-Bruhat’s argument for the wave 
gauge (J6.3D being propagated by the reduced Einstein equations, (16.51) - (16.61) . 
That is, we show that if the wave gauge holds initially, then it is satisfied in 
the whole Cauchy development. To begin, assume the Einstein constraint 
equations and the wave gauge conditions are satisfied by the initial data. In 
order to apply the Bianchi identities in the form Div ff G = 0, we write (16.5|) 
in its equivalent form 

Hij[g\ - i g ij H ap [g\g (Tp = kT.\ j, (6.11) 

and correspondingly the definition of H,j \g\ in (16.71) as 

Hij\il] 2 9ijR(rp[9\9 ^ = [5 1 ] T W ij T dJj,i ~ 9ij 4~ ^p,a) 9 p 

— [ 9 ] 4~ ^i,j 4“ ^j,% 9ij9 ^^(j.p- (6.12) 

Then substituting (16.121) into (16.111) leads to the following equivalent form 
of the reduced Einstein equations 

KTij (-’%j [ 9 ] — ^i,j 4~ kZj,i 9ij9 P^o,p- (6.13) 

Taking the divergence of (16.131) and using the conservation equation (16.61) 
as well as the Bianchi identities in the form Div 9 G = 0, the left hand side 
of (16.131) vanishes and we obtain 

9'i dJj_j r 4“ (f id j t ir 9 ^ id a, pi 4“ Lj dJ a p — 0, 


(6.14) 
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where L? depends only on g and on the Christoffel symbols of g. A direct 
computation shows that the second and third terms in (|6.14j) cancel, so that 
(|6.14j) is equivalent to 

g ap Ui,<r P + L° p uj atP = 0. (6.15) 

The main point now is that (J6.15I) is a homogeneous hyperbolic second-order 
equation for Ui, so that the unique solution of (16.151) vanishes whenever its 
initial data vanishes. Thus, to prove that the wave gauge holds in the Cauchy 
development, assuming the gauge condition (16.31) holds initially, it remains 
only to prove that the derivatives of w* in a direction normal to the Cauchy 
surface vanish initially. For our purposes it suffices to consider the Cauchy 
surface E = {x° = 0} only. We now show that the sought-after vanishing of 
OJi o on E follows from the Einstein constraint equations on E. 

The Einstein constraint equations are given on the surface E by 

Gi 0 [g\=KT i0 . (6.16) 

These are precisely the four Einstein equations for which the second-order 
^-derivatives drop out. Substituting ()6.11|) for the right hand side of (|6.16l) 
gives us 

Gio[g] = H i0 [g} - ^ gioH*p[g]g ap (6.17) 

on E. By (|6.121) . equations (16.171) are equivalent to 

Wi, o + w 0 ,i = giog ap u a)Pl (6.18) 

and contraction with g 13 results in 

g lj (wj, 0 + cuo,i) = ^ 0 g ap ^a, P , (6.19) 

where 6j denotes the Kronecker delta. Now, the wave gauge (16.31) implies 
that uii j = 0 for j ^ 0 everywhere on E, and from this we find that (16.191) 
is equivalent to 

9 ° 3 u<t,q + g 0j uo,o - % g a0 Ua,o = 0. ( 6 . 20 ) 

In the context of approximate locally inertial frames g 00 can always be taken 
to be non-zero, thus the j = 0 component of (16.201) implies = 0 on E. 
This then implies that the j = a component, for a € {1, 2,3}, is given by 

g aa w CT ,o = o, 

and since the metric is non-singular we finally obtain 

OJi, o = 0 (6-21) 

on E for all i € {0,..., 3}. This proves that (|6.21l) holds whenever (16.31) holds 
initially and therefore the wave gauge is propagated by the reduced Einstein 
equations. 
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6.2. The Linearized Einstein Equations in Approximate Locally 
Inertial Frames. In the above subsection, we have shown that the wave 
gauge is propagated by solutions of the Einstein equations (16.5D - d6.6D . as 
long as the wave gauge and the Einstein constraint equations hold initially. 
Based on this, we assume the wave gauge o at the start, and derive 
the linearized Einstein equations in wave gauge in an approximate locally 
inertial frame. This formal procedure accomplishes the following theorem, 
which identifies the non-removable Coriolis terms we seek. 


Theorem 6.1. Assume the background metric c/ij = rjij 4- g^ is given in 
approximate locally inertial coordinates and solves the Einstein equations for 
a source T, while its perturbation g^ = ryj + gij + ehij solves the Einstein 
equations for some perturbed source T = T + eT. In addition, assume that 
g = 0(5), 5 = O(e) and h and derivatives of h are 0(1) as e tends to 0, and 
assume that the perturbation hij is in the wave gauge (J6.3I) . Then, formally, 
substituting the ansatz gtj = rjij + g t] + ehij into the Einstein equations, 
dropping 0(e 2 )-terms and dividing by e, leads to the following linearized 
equations for h: 


K 


-U^hji+Cj^hf-h^RijM = K - -gjig^Trp) --(hj x rf<> - rjjih ap )T ap . 

( 6 - 22 ) 

Here k = — Jj denotes the flat linear wave operator, 

O-r/hji = T] P hji^ p 

and Cjflh) are given by 

Cji (h) = b p lr* pT - h ik S ijkl (6.23) 

where 

Sijki = ^u^kja - r a i a f kj i, (6.24) 

«£. = s i T h - + \ (v^,i +( 6 . 25 ) 


kij — r) (i^ik,j T hjkj hjj k ), (6.26) 

T'kij — 2 i.9ik,j + 9jk,i 9ij,k) j (6.27) 

and we use the notation h lk = g la g kp h C7 p, TO = g lk T k ij and TO = g lk T k ij 
and g l i denotes the inverse of (jij. 

The terms Cjflh), the Coriolis terms, are precisely the terms in (|6.22H 
which do not appear in the linearized Einstein equations in an exact locally 
inertial frame. It is important to note that Cjflh) contains derivatives of 

^Recall that k is negative by our sign convention of the Riemann curvature tensor, as 
we discussed in a footnote in the Preliminaries. Using the opposite sign convention, the 
term —h xk Rijki must be replaced by +h lk Rijki in (16.2211 with no other changes required. 
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gij which are large throughout the approximate locally inertial frame in a 
neighborhood of p, not just at p itself. 

Proof. By assumption, the background metric, g t j = rpj + egg. satisfies the 
Einstein equations for the energy momentum tensor T VJ , that is, 

Rji = Rji[g\ = n (fji - ^gjig ap f a/ ^j , (6.28) 

while the perturbed metric, g t j = gij+gtj+ehij , solves the Einstein equations 
for a likewise perturbed energy and matter source, Tg = Tg + eTij, namely, 

Rji[g + eh] = k ^Tji - ^gjigVPTv^J . (6.29) 

For the purpose of raising and lowering indices, we now introduce approx¬ 
imate expressions for the inverse of g %3 and of gg. To begin with, denote 
the exact inverse of g,g with g and the exact inverse of gij with g . Then, 
setting 

g ij = ri ia v jp g ap , 

we find by cancellation that 

{rj irT -g ia )gaj = {^-g^ivaj + gaj) 

= 5)-g m g*j 

= 5) + 0(e 2 ), 

and since egg is assumed 0(1), we conclude that 

g ij = rf 3 - g ij + 0(e 2 ). (6.30) 

Likewise, defining 

h ij = r 1 ia 'f l 3 P h ap . 

a computation using (16.301) gives 

(r - eh iu ) g aj = - eh ia ) (g aj + eh aj ) = 5) + 0(e 2 ), 

and since gg is assumed 0(1), it follows that 

g ij = g ij - eh ij + 0(e 2 ). (6.31) 

In the argument below we are careful not to take derivatives of (16.311) or 
(16.301) . as this could introduce larger errors in the derivative. 

Now, using the metric ansatz (16.21) and our approximation of the inverse 
metric in (16.311) . we write the Einstein equation (16.291) as 

Rji[g + eh] = k (tji - ^g j ig ap T a p S j - y (. hjig ap - gjih ap ) T ap + 0(e 2 ), 

(6.32) 

which simplifies after substituting Tg = Tg + e Tg and applying (16.301) to 
the form 

= « (Tji - ^g 3 iff p T ap ] + Ke (fji - \rtjirf p f a ^ 


Rji[g + eh] 
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- Y - Vjl h ap ) f ap + 0(e 2 ). (6.33) 

We now expand the left hand side of (16.281) and (|6.33D . Recall that the 
Ricci tensor is the trace of the Riemann tensor, 

Rji = g lk Rijku (6.34) 

and the Riemann tensor is given in terms of the Christoffel symbols by the 
formula 

Rijkl = ^ij[k,l] + g aP (J'ialRpjk ~ ^iak^pjl) , (6.35) 

where we define the commutator 


ij[k,l] 


= r ijk,i - r 




and 


r kij — 2 ^9ik,j 4“ gjk,i gij,k) • 


Moreover, due to (|6.26l) and (16.271) . the Christoffel symbols separate as 


r kij — Tfcij + ef 'kij. (6.36) 

Now, by (I6.34[) and (16.351) . the second-order derivative terms of the Ricci 
tensor are given by g lk and to separate these into its Tjj k - and Tin¬ 
dependence. use that (16.311) agrees with the inverse metric up to 0(e 2 ) errors, 
giving 


g lk ^ij[k,i] — (g lk — eh lk ^j (r im] + ef^*. ^ + 0(e 2 ) 

= g ik ^[k,i] - th ik f ij[kA + eg ik T im] + 0(e 2 ). (6.37) 

To expand the lower-order derivative terms in (16.351) in e, using that all 
terms containing eT squared are of order 0(e 2 ), we obtain 

(T ial F pjk - T iak T pjl ) = {-}i + {-} n + Hm + 0(e 2 ), (6.38) 


where 

{■jt = g ik g ap (r ia iT pjk -r iak T pj i), ( 6 . 39 ) 

{•}n = -e {jf lfp + hd k ff^ {t l(jl f p]k - f l(Tk fp,i) , (6.40) 

{•}lll = eg lk g aP (j'ial^'pjk ~ FicrkTpjl + ^ial^pjk ~ T iak^pjl^j -(6.41) 
Combining (16.371) and (16.38[) . we write the Einstein equations (16.331) as 


Rjl[g d~ eh] — g Rij[ k ,i\ ;] + eg ^ 

+{•}! +{-}n + {-}m + 0(e 2 ) 

= k (fji - ^g j ig ap f ap '\ + Ke [fji - ^ri ] ,ri ap f rjp '\ 


(6.42) 
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To simplify (16.421) further, observe that (16.391) and the first term in (|6.37D 
combine to give the Ricci tensor for the background metric, so that we can 
write the Einstein equations for the background metric ()6.28|) as 

Rjl[g\ = 9 lk ^zj[k,i] + {-}i = K (Tji - ^g j ig ap T (T ^j . (6.43) 

Now subtracting the Einstein equations for the background metric (16.431) 
from the above form of Einstein equations for the perturbed metric (I6.42p . 
we obtain 

e 9 lk ^ij\k,l] ~ eh lk t ij[ k: i\ + {-}n + + 0(e 2 ) 

= ^ (f jt - h, 3l !fPf\ - ^ ( h]l JfP - g jlh "P) f ap . (6.44) 

For further simplification of (16.441) . we use the formula (16.351) for the 
Riemann tensor of the background metric, which leads to 

- + {-}n = —eh lk (r y[ fcj q + g ap {Tiai^pjk ~ ^iak^pji) ) 

- eg ik h° p (Y ial T pjk -T i(jk T pjl ) 

= -eh ik R ljkl - eg ik h ap (T ial T pjk - T lak T pjl ) 

= -eh ik R ijk i-eh ik S ijkl + 0(e 2 ), (6.45) 

where we substituted (16.301) for g lk in the last equality together with the 
definition of Sij k i, (I6.24p . Thus, substituting (|6.45l) into (|6.44l) . the linearized 
Einstein equations obtain the form 

eg ik t ij[kA + {•}„! - eh ik R ijkl - eh ik S ijkl + 0(e 2 ) 

= Ace (f jt - - y ~ 9jih ap ) f ap . (6.46) 

To derive the final form of the linearized Einstein equations in approximate 
locally inertial coordinates it remains only to impose the wave gauge condi¬ 
tion. The main step is the following lemma. 

Lemma 6.2. Assume hij satisfies the wave gauge (16.31) . then 

eg ik t im] + {}m = + eb p l TV + 0(e 2 ), (6.47) 

where b p J a is defined in (16.251) . 

Before giving the proof of the Lemma, note that (16.471) implies the lin¬ 
earized Einstein equations (16.221) . by substituting (16.471) into (I6.46P and di¬ 
viding the resulting equation by e, yielding 

+ b p l f- h ik S ijkl - h ik R ljkl + 0(e) 

= K (Tji - f (hurfP - r] jl h ap ) T ap , (6.48) 


which is (16.221) . To complete the proof Theorem 16.11 it remains only to prove 
the lemma. 
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Proof of the Lemma: We first show that using the wave gauge condition, 
(HCU), the second-order derivative terms combine to the flat wave operator 
□, 7 . To begin with, substitute (|6.30D for g lk in the second-order terms in 
(|6.44l) and use that g V] = 0(5), we then get 

eg ik r ij[k ,i] = e(rt k --j k )t im] + 0(e 2 ) 

= erf k Y im] + 0(e 2 ), (6.49) 

while a straight forward computation leads to (c.f. (16.81) 1 

r.ikf' _ f'k 

V ij[k,l] jk,l jl,k 

= 2 — h k j,ki ~ h k ij k + r] ap hji )a p'j 

= tpnhji - \cH (h\ k - \h k k \ - Idj (h\ k - \h k k A (6.50) 

where we raise indices on hij and T,j k with the Minkowski metric. From the 
definition of Y t j k , (I6.26|) . a straight forward computation shows that 

h\k-\h\i = rf p ti. P , (6-51) 

and substituting f|6.51[) into (I6.50P yields 

ri ik f~ \d t (rf p tj* P ) ~ \dj ^f p f iffp ) ■ (6.52) 

Let us remark at this point, if locally inertial frames existed, we could 
impose the gauge condition ri ap Ti ap = 0 , which would be propagated by 
the resulting flat wave equation, in which case the right hand side of (|6.52j) 
would reduce to However, in an approximate locally inertial frame, 

we work with the wave gauge condition (16.31) which is propagated by the 
full non-linear Einstein equation. In order to use (16.31) in (16.521) . we have 
to separate r] ap Ti ap from the remaining terms in the wave gauge condition. 
Naively, this could be achieved using the approximate inverse of gij in (I6.3ip . 
then the wave gauge condition assumes the form 

V ap r iap = (g° p + eh ap ) t iap + 0(e 2 ). (6.53) 

However, using (16.531) in (16.521) . at this stage we would loose control of the 
error due to the derivatives d[ and dj on the gauge functions. Therefore, we 
approximate g 13 in a way such that we can control derivatives of the error 
terms. For this, we introduce 

eh 13 = g 13 — g 13 , (6.54) 

g ij =g ij ~rf 3 , (6.55) 

from which one obtains the exact identity 

g ij = g ij + eh ij = r] ij + g ij + eh ij . 


(6.56) 
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Now, eh 13 and g 13 are both of order e because substituting (16.311) into (16.541) 
gives 

eh ij = -eh ij + 0(e 2 ) = O(e), (6.57) 

and similarly approximating g v in (16.551) by (16.301) leads to 

~g ij = -g ij + 0(e 2 ) = 0(e). 

Regarding derivatives of g lJ , there is no reason to assume that g p) j = g l \ is 
0(e), so we must keep all such terms in the linearized equations. However, 
eh lJ t = O(e), since the derivative of the inverse metric is given by 

g 3 g = -9 ia g jp 9a P ,i, 

which implies by (16.541) that 

eh ij tl = <h {g' j ■ - f /') 

= -9 ia 9a P ,i9 jp + 9 ia 9a P ,i9 jp 

= 9 m 9 3P (9a P ,i ~ 9a P ,i) ~ e(g ia h 3P + h lr7 g 3p ^j g ap%h (6.58) 
and using now (16.571) and 

9a P ,l 9cr P ,l £h a p i 0(e), 

we hnally obtain 

eJA, = 0(e). (6.59) 

From this we conclude, that the inverse of g^ in the form (16.561) separates 
off the 0(1) term if 3 from the 0(e) terms eh 13 and g lJ , and (16.58[) confirms 
that raising the indicies on h keeps h and derivatives of h order 1. 

Using (I6.56p . we write the wave gauge condition in its equivalent form 


v ap r i(Tp = -(r p + eh^)f 


tap- 


(6.60) 


Now, substituting (16.601) into (16.52[) we obtain 


l cr P J i 
(6.61) 


V ik r im] = -f^hn + -di (( g ap + eh"*)Y jap ) + -dj ((r p + eh^Y, 

so that using (16.591) . g‘ J t = g t] { and eh 13 = O(e) = g 13 finally leads to 

eri ik r ij[kfl = + e -g ap t jap + f, i* P + 0(e 2 ). (6.62) 

To finish the proof, we combine the derivatives of g %3 in (16.621) with {-jin 
and substitute (16.301) for g lk which yields up to 0(e 2 ) errors that 


{.} m+ e -g" P l t Jap + e -r 3 T lap 


= e r* ,r 


i 


al L jk 


p/c pa I p/c per p/c per . ^ *-apy\ , - 

^ akjl ' ^ aljk ^ ak ^ jl ' r\9 jap \ 


al L jk 
rr-n/c 


crpp 


l ap 


1 


5JY p a - 51 f p + 5 P t T al - 5 p 6JY k k + - + r^sTA ^ r , 


= e 
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where we raise indices on T and f with the Minkowski metric. Comparing 
the terms in the braces of the previous equation with the definition of the 
coefficients in (16.251) immediately yields 

{'}m + (6-63) 

Finally, combining (I6.63|) with (I6.62|) . we obtain (16.471) . This proves the 
Lemma and completes the proof of Theorem 16.11 □ 

6.3. Coriolis Accelerations for Gravitational Waves in Approxi¬ 
mate Locally Inertial Frames. Our goal here is to prove that in each 
approximate locally inertial frame there exists gravity waves in the wave 
gauge which solve (16.221) . such that Cji(h ) is of order 0(1) in a neighborhood 
of p. In this case, Cji(h ) represents an acceleration to the gravity wave 
which we identify as a non-removable GR Coriolis acceleration in analogy 
to the classical Coriolis force. 

We see from (16.221) that even though g is only Lipschitz continuous and 
the linearized Einstein equations contain second derivatives of g , the delta 
functions cancel out in every approximate locally inertial frame because they 
only appear through the curvature tensor. Due to cancellation, the curva¬ 
ture tensor does not distinguish an approximate locally inertial frame from 
an actual one. This then begs the question as to whether it is possible for 
the first-order derivatives to similarly cancel in approximate locally inertial 
frames. Asked differently, do the first-order derivatives of g distinguish ap¬ 
proximate locally inertial frames from actual ones? The purpose of this sec¬ 
tion is to answer this question in the affirmative by isolating a non-vanishing 
scattering effect in approximate locally inertial coordinate frames which dis¬ 
tinguish them from actual locally inertial frames, and thereby complete the 
proof of Theorem 15.11 of the Introduction. 

Theorem 6.3. Under the assumptions of Theorem [6771 in each approximate 
locally inertial coordinate system at p, in the sense of Definition \3.1\ there 
exists indices j, l, a, p, t with p / t, such that 



Moreover, within the class of initial data satisfying the wave gauge 
there exist solutions h of (16.221) and l,j from, (|6.64D such that 

\Cji(h)l > hi, (6.65) 

in some open subset of U. 

The theorem makes things simpler than one might expect because it ap¬ 
plies whenever any derivative of g fails to vanish. 

Proof. The well-posedness for (16.221) follows from the standard existence 
theory for linear hyperbolic PDE’s of second-order m e]. It remains to 
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prove (|6.64j) and to show that there exists initial data for which the resulting 
solution satisfies (16.651) . 

We first prove (16.641) . By (13.31) of Definition 13.11 there exist indices a, /?,7 
for which \g a p,-y\ > M in U. Now, since g a p,i can be expressed as 

9aj3 ,7 — bo 7 7 -f- Tfa a , 


we know that 

| r'cK/3'y 4“ b/5 7 n: — .4a:/1.7 O Af, 

yielding that there exists some set of indices, which for simplicity we again 
label as a, (5, 7 , such that 


\*%\ > M/2. 


( 6 . 66 ) 


We next verify that (16.661) alone implies (|6.64j) . To start, recall the defi¬ 
nition of the coefficients l /‘/ a , (16.251) : 


yi = SJt p ia - + 8?th - 8^t k ak + 


■ J* 7 


' jl 


~pT . ~pr 

Vj*9 j + T]ia9 J 


To prove that (16.661) implies (I6.64p . we compute b r "j a under four different 
conditions on the indices l,j,a,p,r. For these explicit cases, we do not 
use the summation convention, but rather calculate bj/ a for fixed indices 
l,j,a,p,r when explicit relations between them are assumed. 

Case I: Assume er = t, p / a, a £ {j, l}, and p £ {j, l}. Then by (|6.25D 


m = 0 - r£ + 0 - 0 + io = t p jV (6.67) 

Case II: Assume p = l, t = a, j ^ l, and o £ {j, l}. Then by (16.251) 

b l jia = 0 “ d°T' l j l + 0 — 0 + -0 = —T l ji- ( 6 . 68 ) 


The constraints on the indices in Cases I and II together with the symmetry 
ff. = rj allow us to solve for any D) on the right hand side of either (16.671) 
or (I 6 . 68 p . except for the case when p = j = l. To address this last possibility 
we require one additional case: 

Case III: Assume p = l,r = a,a^l=j. Then by (I6.25P 


b% = 0 - £f J, + 8\Yh - 0 + io = -t l u + TZf (6.69) 

Now assume rjg > M/2 for some values of a,/3 ,7 such that not all of 

cc,/3, 7 are equal, then, (since T is symmetric in the lower two indices), 
either a ^ {/3, 7 } or /3 ^ 7 , a = 7 . If the possibility a ^ {/3, 7 } holds, then 
we can apply (16.671) of Case I to conclude 


I b 


acr I 

£ 70-1 
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for some u/a, (again, we do not sum over a in the above equation). If the 
possibility /3 ^ 7 and a = 7 holds, then we can apply (16.68[) of Case II to 
conclude 

\KU = |f£J > \m 

for a {a,/3}, and hence for some a ^ a. This confirms 16.641 for every 
possibility except the case a = /3 = 7 . So assume a = f3 = 7 . Using this in 
(16.691) of Case III with k / a, we conclude 

Ka = rL-K k ak - (6-70) 

Now, substituting (16.681) for some cr ^ {a, k} and k 7 ^ a, (that is, b^ k(J = 
—f^ fc ), into (16.701) yields 

fSa = -&-&. (6-71) 

which immediately implies 

fcl + K k ak \ > fc + b « k ak | = |f« a | > y. (6.72) 

From this we conclude that either | > M/4 or &l > M/4. In 
Summary, we have proven that (16.641) holds for some coefficient bjj a such 
that p / t. 

To verify (16.651) . we need to prove that for some j, l, there is not sufficient 
cancellation among the expression for Cji given in (16.231) to make \Cj/\ < M/4 
for all h. To prove this, recall that by (16.641) . there exists fixed indices 
j, l, a, p,T (with p 7 ^ t) such that 



Now, we first decompose Cji as 




(6.73) 


where we define Q to be the sum over all terms on the right hand side of 
(16.231) . which do not involve T/ t , for the fixed indices cr, p , r, divided by the 
non-zero b^ a from (16.641) . 

Now, to finish the proof, fix an arbitrary point q £ U. Then, by (|6.73p . 
it remains only to show that there exist an h such that, at the point q, T^ t 
satisfies 


M 
8 b 


Q < Tpr < 


M 

8b 


Q, 


(6.74) 


for b = \bjj a \. Indeed, (16.741) is the sought after estimate (16.651) at the 
point q £ U, which by continuity would hold on some neighborhood of 
q. Therefore, it remains only to show that there exists enough freedom to 
choose h and derivatives of h at the single point q to make Y a satisfy (16.741) 
at q. For this final step, the only restriction among all T\- are the Einstein 
constraint equations and the wave gauge condition (16.31) . because fT at q 
can be considered to be initial data from which h is determined. 
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The rigorous analysis of the constraint equations is beyond the scope of 
this paper, but all we are assuming here is that assigning some value for 
the T£ t , at the single point q, is consistent with the well-posedness of the 
constraint equations. 

Finally, by |2], the solution hij is determined by initial data h a p and 
h a ^o for a, j3 € {1, 2,3} on £ = {x° = 0}, leaving the freedom to arbitrarily 
assign hto and /ijo,o for i € {0, ..,3}. These are determined by the wave 
gauge conditions (16.31) in its equivalent form 

3 

^9 {hikfl T foo,fc ' 9 

a,/3= 1 

where the right hand side only depends on spatial components of h. Thus, 
the wave gauge conditions (16.31) can be satisfied by choosing /ijo,o accordingly 
(c.f. p. 164, chapter 7, in [2]), which then leaves /i Qj g and h a pfi for a, (3 € 
{1,2,3} free to assign as initial data and this freedom clearly suffices to 
arrange for (I6.65P to hold. In summary, we have shown that one can choose 
the value of Y a pT at q to satisfy (|6.74l) at the start and then it is consistent 
to solve the wave gauge and the constraint equations simultaneously with 
that value. This completes the proof of the Theorem. □ 


7. Conclusion 

The point of departure for this paper is the authors’ recent work in mm 
which gives the first proof that spacetime is locally inertial at a point of 
shock wave interaction. The method in mm is tailored to the simplest 
case of shock wave interaction, to the case of interaction between two shock 
waves from different characteristic families in spherical symmetry, and it is 
an open problem whether regularity singularities exist for more complicated 
shock wave interactions. We know of no physical principle that can rule out 
regularity singularities ahead of time, and only a mathematical proof can 
resolve the problem as to whether regularity singularities exist. In this pa¬ 
per we clarify and motivate the open problem by investigating the physical 
implications of regularity singularities should they in fact exist. To clarify 
the issue, we make the distinction between the smooth atlas of C 2,1 coordi¬ 
nate transformations and the larger atlas of C 1,1 transformations. We then 
prove that, restricting to the C 2,1 atlas for a metric that is Lipschitz con¬ 
tinuous at a point of shock wave interaction, the closest one can get to a 
locally inertial coordinate system is one which is approximate locally inertial 
in a sense we clarify. We then linearize the Einstein equations in an ap¬ 
proximate locally inertial coordinate system, and identify and characterize 
the Coriolis type terms which we prove will only vanish in a true locally 
inertial coordinate system, should one exist. The open problem of regularity 
singularities, then, is the problem as to whether the approximate locally 
inertial coordinate systems can be improved to locally inertial coordinate 
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systems within the larger C 1,1 atlas. If locally inertial coordinates do not 
exist within the C 1,1 atlas, then the scattering of gravitational radiation by 
a regularity singularity would produce quantifiable physical effects analo¬ 
gous to non-removable Coriolis type forces, and these Coriolis effects are 
characterized in this paper. 
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